function e = expected_cost(chat,cbar)

% hi = @(z) z.*cbar.*exp(-cbar.*z) ;
% 
% for i = 1:size(chat,1)
%     for j = 1:size(chat,2)
%         e(i,j) = integral(hi,0,max(chat(i,j),0))./(1-exp(-cbar*max(chat(i,j),0))) ;   
%        if (1-exp(-cbar*max(chat(i,j),0)))==0
%            e(i,j) = 0 ;
%        end
%     end
% end

k = max(chat,0) ;
e = ((-k.*cbar.*exp(-cbar.*k)-exp(-cbar.*k)+1)./cbar)./(1-exp(-cbar*max(chat,0))) ;
e(isnan(e)==1)=0 ;

% 
% d = cbar ;
% c = k ;
% 
% exp(-c.*d).*(exp(c.*d)-c.*d-1)./d